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RADIATI VELY DRIVEN WINDS FROM MAGNETIC, 
FAST- ROTATING STARS 

BY 


Steven Nerney 


ABSTRACT 


An analytical procedure is developed to solve the 
magnetohydrodynamic equations for the stellar wind 
problem in the strong-magnetic field, optically thick 
limit for hot stars. The slow-mode, Alfven, and fast- 
mode critical points are modified by the radiation terms 
in the force equation but in a manner that can be treated 
relatively easily. Once the velocities at the critical 
points and the distances to the points are known, the 
streamline constants are determined in a straight-forward 
manner. This allows the structure of the wind to be 
elucidated without recourse to complicated computational 
schemes . 
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I. 


INTRODUCTION 


Since the study of Weber and Davis (1967, WD) showed 
that the sun could be spun down by the solar wind over 
its main sequence lifetime, work has continued on angular 
momentum loss in the solar wind and in stellar winds 
(Pizzo, et al., 1983; Nerney and Suess, 1975; Nerney, 
1980; Goldreich and Julian, 1970; and Friend and 
MacGregor, 1984). 

Early-type stars present additional problems in 
understanding angular momentum loss as stellar winds from 
0 and B stars are primarily driven by their intense flux 
of ultraviolet radiation. The radiation has an 
associated outward force due to Thomson scattering by 
photospheric electrons as well as a line radiation force 
that acts on metals in the wind. Unfortunately, 
analytical solutions are not possible as the basic 
stream-line constants are not known until a full 
numerical solution is generated. The complexity arises 
because the force due to line radiation depends 
nonlinearly on the acceleration in the wind. The 
critical point analysis shows that there is a focus of 
singular points which the numerical solution must be 
tangent to. This is accomplished through an iteration 
scheme that produces an acceptable physical solution. 

The purpose of this study is to attempt to 
analytically model radiatively driven stellar winds with 
strong magnetic fields from fast-rotating stars in the 
opt ically-thick limit. This will allow a simplification 
of the problem so that it becomes readily apparent how 
the streamline constants determine the solution. Also, 
the modified magnetohydrodynamic critical points are 
treated and a procedure is developed for constructing 
wind models which does not require a complicated 
interaction scheme. 
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II. EQUATIONS OF MOTION 


The steady, spher ically-symmetric, dissipationless, 
magnetohydrodynamic (MHD) equations, together with terms 
that give rise to radiative bulk forces may be written as 
follows : 

2 # 

(1) 4ir pv = M, which is the continuity equation 
where M is ?he mass-loss rate in kg/sec, r is radial 
distance from the center of the star in the usual r,0,<f> 
spherical polar coordinate system, p is the mass density 
in kg/m , and v r is the radial velocity of the wind. 


(2) 



1 dP 
p dr 


GM (1-r) . 

2 

r 





8ir p r 


_ d_ 
2 dr 


(rB.) 2 + 
<P 


which is the radial momentum equation where P is the 
pressure, G is the gravitational constant, M is the mass 
of the star, r is a constant (see below) related to the 
radiative force due to Thomson scattering of continuum 
photons by electrons, v is the rotational (azumuthal) 
velocity of the wind, B^ is the azimuthal component of 
the magnetic field, and <p f L is the force/kg due to line 
radiation. The term in B is the radial component of 
the j x ^ force. ^ 


Q0 Hi 

(3) r = -r-- -— — where L is the stellar luminosity, c is 

the speed of light, and P e is the electron scattering 
opacity per unit mass. 


(4) 


f. = 


tgm 


kt 


-a 


which is an approximate law derived in the original 
theory of Castor, Abbott, and Klein (1975, CAK) where k 
is a constant which is a measure of the mixture of the 
number of strong radiation lines in the wind, t is an 
optical depth parameter, and oi is a constant which 
depends on the importance of optically thick and thin 
radiation lines and is model-dependent. 
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(5) t = o g p v ^ (^—) and v th the thermal velocity 

of the ions which absorb and scatter radiation. CAK used 
values of k = 1/30 and a = .7, but I will use a = 1 and 
choose various values for the constant that k appears 
in. This Sobolev high-velocity gradient approximation to 
the line radiation force is discussed in Mihalas, 1978, 
p. 561 . 

The azimuthal equation of motion and the induction 
equation (a form of Faraday's law) are the same as in 
Weber and Davis (1967, WD) and Friend and MacGregor 
(1984), namely: 


(6) 37 (r V = 3777 - 37 (r V 


(7) <v x E) t = 0 = i [r („ r B + - v^B r )] 


as V • S = 0 gives 


2 

(8) B r ~ f t/ r f° r spherical symmetry and as B r /P v r is a 
constant, (b) integrates to 1 


( 9 ) 


rv , - 

<p 4npv 


rB, = L where L is the stream constant 
<P 


total angular momentum per unit mass which is the sum of 
the angular momentum carried in particles, rv , as well 
as that in the electromagnetic field. Faraday's law 
integrates to 


(10) r(v r B, - v.B ) = constant = - 2 f fa and (9) and (10) 
together can be solved for v, and B, as in Weber and 
Davis (1967). * * 


V — V 

Qr rA r 

(11) % = t. ~ 


rA 1 - M 


A r 


XXXli-3 



( 1 2 ) B a 
9 


B 


Qr 


2 2 
r A ' r 


r VrA r A 2 (1 - "aV 


where is the rotation rate of the region where the 

field lines are rooted (usually taken to be the 

photospheric value), f. is the conserved radial magnetic 

flux, v A is the radial velocity at the Alfven 

radius ( r - ) where v r equals the speed of an Alfven wave 

in the medium, and M, is the radial Alfven mach number. 

Ar 


(13) M a 2 
A r 



p 3 /p . 


Closure is obtained by using the polytropic 
approximation. 


(14) P = P (p/p ) y where P and p Q refer to values at 
the base of°the w?nd (whether the base is in a corona or 
photosphere), and y is the polytropic index. 

This approximation is often used when the heating 
sources in the wind are not well-known. Generally, 
modellers replace this with an isothermal approximation 
for hot stars and use the stellar effective temperature. 
This is formally incorrect as it leads to inappropriate 
asymptotic states of the wind, Hundhausen (1972, p.9). 
However, the use of effective temperatures of even 50,000 
K as for 0 stars does not lead to significant problems as 
it would for stars with coronae. Most early-type stars 
show x-ray emission. Therefore, I will examine a hot- 
wind solution and the polytropic approximation is 
appropriate . 

Th^ radial momentum equation, (2), can be simplified 
using a = yp/p (a is the sound speed) together with 


(15) f T = — — =r • Using a = 1, (2) eventually 

, J_j , 01“ 1 Q IT 

reduces to p 
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dv 


(16) J- 

r 


r 

dr 


, 2 . 2 . r 0 2 . 2 
< Y r ~ A r 1 t 2a * % ~ 


< 1 - r >) + 2 v r % A r \ 


, 2 . 2 * , 2 « , , » 2 , , 2.2 
( v - A ) [v (1-A) - a J - v A 
r r r r 9 


where A 2 = B 2 /4ttp, A 2 = B 2 /4rrp, 
r r <j> <j> 

2 2 * GM 

azimuthal Alfven speeds, v g = — 
at r, and r 


the radial and 


, the escape velocity 


(17) A = (— 


4 7T 


a 


M 0 v., 

e th 


■) r GMk or, for a = 1, 


(18) A = 


kL 


CM v 


th 


A is the essential radiation parameter in this wind 
model and can be estimated for a particular star. For 
instance. Friend and MacGregor (1984) model the 06ef 
:ephei with a = .7, M » 5.2 x 10 - 6 M o /year , 


star X Cephe: 


T = 6.5xlO^K, L = 6.76xl0^L o , and k was specified by 
using the CAR value of 1/30. These numbers give A = .66 
in the current model. It is difficult to estimate A 
because k is model dependent, v t ^ depends on both the 
base temperature as well as the mass of the ions that 
absorb and scatter radiation, and mass-loss rates are not 
well-known for many stars. 


Therefore A will be specified for a range of values 
and a variety of models will be calculated. 


The general results of this study are based on the 
strong magnetic field limit of the MHD equations. This 
allows the explicit calculation of the position of the 
critical points based upon a similar prescription for 
non-radiat ive strong-field winds in the study of Hartmann 
and MacGregor (1982, HM) . 
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III. THE CRITICAL POINTS 


A . Critical V e lo cities 


The WD theory gives rise to three critical points in 
the radial momentum equation. These points correspond to 
distances from the star where the speed of the stellar 
wind equals the speed of the characteristic MHD 
disturbances that occur in the fluid approximation; 
namely, the slow-mode wave, Alfven wave, and fast-mode 
wave. Setting A = r = 0 for the moment, the denominator 
of equation (16) goes to zero for three different wind 
speeds and, hence, the numerator must be zero at these 
critical points. These boundary conditions partially 
specify the physical characteristics of the wind. In the 
current study, radiation changes the nature of the 
critical points in a manner that can be explicitly shown. 

2 2 

The strong-field limit requires that v^ << A^ at 
the modified slow-mode point, r s . This is equivalent 
to m A r « 1 at r g . 


This will be true for sufficiently strong magnetic 
fields and is justified a posteriori. Getting the 
denominator of (16) equal to zero: 


A 


r 


2 



(1 - A) 



= 0 


Faraday's law 


can be written as 


v 



2 

4 _ 


2 

(v - nr) 

<p 


and as v^ is equal to the corotation value at r g to 

2 2 2 2 
0(M Ar^ we can neglect v r A /A^ in the strong-field 

limit. Therefore, ^ 


(19) 


rs 


1 -A 
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which reduces to a well-known result when 


At the Alfven point, the radial velocity is assumed to be 
much greater than its value at r g so that 


v 


rA 


(1 - A) » a‘ 


Setting the denominator of (16) equal to zero. 

A 2 

(20) v = A + at r at which reduces to the well- 

k r IT 1 A cl 

nown 

result v^ = A at r A when A = 0. 

The mod^f ied 2 fas^-mode velocity is found by 
assumming A (M Ar > ^ at r f* In a s i- m il ar manner, 



B . Optically-Thick Slow-Mode Point 
Expanding Equation (11) for v^ 

2 

(22) v = flr(l - M Ar + £— 2 M Ar + . 

r A 


) 


For clarity, let v, = nr(l - 5) where 6 < 1. 
Setting the numerator equal to zero at r s 

2a 2 - v e 2 /2 + (v^ - 2V c A + /A r ) = 0 . 

Using equation (10) gives 
2 

(23) 2a 2 (1-r) + O 2 r 2 = 0 to 0(6 2 ) 

as terms of 0(6) exactly cancel. 
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Using the non-dimensional parameters in HM, 

2 2 3 

Z = r /r , ft = a GM/r , and the non-dimensional 
Parkercr it ical point Z = ~ , we find 

P Or- a A 


(24) Z s ~4 + “ 2 2 p z s 3 - 2 p (1 - r) - 0 
a 

o 

which reduces to equation (9) in HM for the isothermal 
limit. Equation (24) does not explicity depend on the 
radiation parameter, A. The distance to the slow-mode 
point is the same in the radiative wind (except for the 
hottest stars where r is not negligible) as in the non- 
radiative MHD wind, except that the radial velocity in 
the radiative case is greater by the factor 1/ 1 - A. 

Another point is that the calculation of Z assumes 
1 s 

v^ is given by the corolation value at r g . This is 
accurate to 0(M ), yet due to a fortuitous cancellation, 

AIT 2 

the equation for Z g is accurate to 0(M Aj _) so that these 
results and HM's are more accurate than anticipated. 

2 2 

To solve (24) for Z s , we must know a /a which 

2 ° 2 

depends on v rQ . In the isothermal limit, a /a Q = 1 
and this problem does not arise. However, 


(25) 




(26) B 


2 (y - 1) 

Y + 1 


A second equation in Z and v can be found by 


ro 


examining the energy/kg in the flow. 


(27) E = j (1 


,, 2,1 2 

A) v^ + - v , + 


a 2 v e (1 ~ r, _ nr R r % 
2 '<t> ' Y-l 2 V 
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The Poynting flux term can be calculated from 
equation (9) so that 


, 2 v * 

(28) E = -i- (1 - A) v + -j v^ + y~ + QL - ar 


Consider the two terms in v. using v, = £2r(l - 6) 

<P <P 


1 2 

2 % " flr % = 


.2 2 0 
- fl - s r - (1 - 5 2 ) 


so that terms of 0(M Ac ) have again cancelled. 

Settinq E = E we find after some work 
a o s 


(29) (l^A) - (1 + -ij.) 


a z 
o s 


2z p ((1 - f) (1-r) -§- (2 S 2 - D) 


This reduces to equation (13) in HM with y = 1 and A = 0 

a 2 2 

provided we replace j- by the isothermal term a £np. 


Equations (24) and (29) can be self-consistently 
iterated between to find v and Z g using HM's initial 
guess, modified for radiation 



3Z p (1-r) 

1+3 a 2/3 Z (l-r) 2/3 


. 2 2 

which assumes a /ao » 1. This value for Z allows the 
calculation of v in (29). The iteration procedure 
rapidly converges to a solution. 
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C . Optically Thick Fast-Mode Point 


We have shown that tor M^ >> i v 

a.L rt 

which together with equation (12) gives 


_ 


l-A 


v 


rf 




3/2 2 

v ' r. 

m (1 - ^ 


v 


TJ2 


rf 


where v m is the usual Michel velocity, Michel (1969) 


2 2 

(30) v 3 = so that 

m • 

m 


(31) v 


m 


r£ ( 1- A ) 1/3 


In order to calcuate z^, we need to know v f which 
can be calcuated from (11): ^ 


v 


v f — { 1 — e } « — (1 - — — } » e < 1 
<))f r 1 J r ^ v r f 


which is formally correct with or without radiation. 

Setting the numerator of (16) equal to zero at Z^: 


(32) 2a 2 - v 2 /2 (1-r) + v. (v. + 2 (v - L/r)) = 0 

e tp <p (|) 


u — i 

The term in v. reduces to — ^ (1 - 45) to 0(M. ). to 


A 

Equation 9 ( 32) now redn6es to 


Ar 
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2 2 , 

a f 2 2 4 2 v m ^/^ 

(33) Z* - Z Z. (1-D + a Z^ Z a - 2Z_ -Ar (1-A) = 0 

2f p f P A A 2 


where 


(34) 


v ro (1 - A) 

Z 2 v 
r m 


1/3 


Equation (33) reduces to (17) in HM in the 


isothermal limit with A = 0 except for the term in — j . 
This is 0( M Ar ) and should not appear as other term© of 
the same size have been neglected in calculating Z^. 

As Z A appears in the calculation of Z f we must 
generate another equation for Z^ and Z^. The Alfven 
radius is not important except insofar as it is necessary 
in the calculation of Z f . If the wind passes through Z s 
and Z f , it must of necessity pass through Z A , Goldreich 
and Julian (1970). 

Settinq E = 

o f 


(35) 


m 


2a 2 Z 


2 , 2 , 
a f /a o _1 
Y-l 


v 


ro 


2a 


with 


+ h (1 - !-) (l-D + 7 


(36) n i | (1-A) 1/3 + 1 
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A good first approximation to Z A neglects l/Zf and 
2 2 

sets a^ /a » 0. This allows the calculation of Zg in 
(33) and the iteration proceeds as usual. 
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Ill . 


CONCLUSIONS 


While the results of this study are preliminary, 
several important conclusions can be drawn from the 
analytical techniques that have been developed. The 
strong-magnetic field, optically thick limit of the MHD 
equations produces physical solutions through a tractable 
technique that shows both how the critical points are to 
be treated and gives the values of the velocities at 
these points as well as the distances to the critical 
points . 

The streamline constants are determined once the 
solution has been generated at the critical points. In 
particular, the distance to the slow-mode point and the 
radial velocity at the base must be iterated on to 
produce physical solutions. Once this is done, the mass- 
loss rate can be determined from the value of the density 
in the photosphere together with the radius of the star. 
Additionally, the total energy/kg is also determined. 

The angular momentum/kg is determined after the iteration 
procedure to determine the Alfven radius and fast-mode 
radius is completed. 

Preliminary numerical solutions have been generated 
and are being compared to the work of Friend and 
MacGregor, 1984, although it is still too early to report 
these results. These solutions have been found by using 
an IBM-PC and, in principle, could have been done on a 
hand calculator. 
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